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We demonstrate our simple strategy for renormalization with QED at one-loop level, basing on 
an elaboration of the effective field theory philosophy. No artificial regularization or deformation 
of the original theory is introduced here and hence no manipulation of infinities, ambiguities arise 
instead of infinities. Ward identities first come to reduce the number of ambiguities, the residual 
ones could in principle be removed by imposing physical boundary conditions. Renormalization 
Q ' group equations arise as "decoupling theorems" in the underlying theory perspective. In addition, 

, a technical theorem concerning routing of external momenta is also presented and illustrated with 

■ the self-energy and vertex function as examples. 
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INTRODUCTION 



All the presently known quantum field theories for practical use are beset with ultraviolet divergences 
to certain extent in perturbative frameworks. These divergences imply that such theories are not well 
defined at short distances before renormalization is carried out, at least withinperturbative frameworks. 
The celebrated BPHZ program is the standard algorithm for renormalization[l|-|3[ that is widely employed 
by particle and field theorists, which can lead to a mathematically well-defined perturbative formulation 
of a field theory in the end. But the procedures involved are not easy to be endowed with a natural logic 
that would be readily to accept. In this connection, Epstein and Glaser (EG) had already shown that 
finite perturbation theory of quantum field theories can be obtained through constructive methods basing on 
causality conditions [Hill where Feynman rules are only valid at tree level, which provides more mathematical 
rationales for conventional procedures like BPHZ [g, Hi- Recently, other mathematical structures hidden in 
the perturbative formulation (especially in BPHZ) of quantum field theories has become a rapidly growing 
area of research (8rll0|. which implies that there are more regularities worth exploring under the veil of 
renormalization. 

From more physical point of view, according to the effective theory philosophy, the appearance of ultraviolet 
divergences can be understood from the fact that the presently known formulation of a quantum field theory 
is necessarily an effective one with the sophisticated structures or parameters dominant at short distances 
or larger energy scales being ignored. Unfortunately, the sophisticated structures underlying the present 
formulation of quantum field theories are actually unavailable or not discovered yet. Nor could us simply 
guess about the exact answers. The divergences or the ill-defined objects are the prices we paid for such 
forced ignorance or simplification. This is the start of the effective field theory philosophy that is now widely 
accepted and applied by physicists. Fortunately, this effective field theory view also implies that the ignored 
short-distance details could at best affect the physical behaviors of effective field theories through local 
interactions or operators that could be parametrized as corrections to the actions of quantum field theories, 
which, after being absorbed into the normalization factors for effective field theory's couplings and operators, 
leads to the renormalization of quantum (effective) field theories. Otherwise they would be "captured" as 
explicit and active degrees of freedom. 

In practical computations, one usually has to first introduce various regularization schemes and interme- 
diate divergences to be removed later, which are supposed to imitate the roles played by the true underlying 
structures or parameters. Actually, with such operations, the perturbative expansion only takes rather for- 
mal meaning due to the presence of infinite counter-terms. Moreover, for the regularization schemes or 
deformations to be legitimate and useful, they must not harm the bulk structures of the original theory as 
well as being self-consistent. In particular, important symmetries and global properties like the topologies 
of spacetime and Hilbert space, unitarity and causality, etc., which serve as foundations for a quantum field 
theory, should be preserved as much as possible, at least no drastic or uncontrollable changes should be 
introduced. However, the status of the principles just stated are extremely hard to precisely or profoundly 
assess within any known regularization. As a matter of fact, the comprehensive and global structures of 



2 



a well-defined quantum field theory have never been fully delineated, and we heavily rely on rather formal 
evaluations to deal with this very difficult issue. Sometimes, to pursue efficiency, even the logical foundations 
of the theory would be severely mutilated [llj]. Thus, we only accumulate our faiths in one regularization 
method through intensive tests against experiment data. Up till now, there exists no regularization that 
could be satisfactorily in service in all contexts, each has certain shortcomings that will fail somewhere, or 
alternatively, no one could prove to work everywhere 1 . For example, the celebrated dimensional regulariza- 
tion works excellently for standard model, especially for the strong interaction, but not for supersymmetric 
theories at all jl2Hl9| . Of course, other schemes also have their own shortcomings, for a comprehensive 
discussions of the advantages and disadvantages of various existing regularization methods, please refer to 
Ref. ■ This is also an indirect reflection of our ignorance of the grand structures or contents of quantum 
field theories in depth. Therefore, most regularization schemes are actually used in a superficial or formal 
sense, only appearing in the intermediate stages of calculations. After the artificial and transient excursion to 
parametrize our ignorance of the underlying short-distance physics, we "return" to our "original" quantum 
field theories with the deformations removed somehow. This practice is quite established within perturbative 
regimes for many physically interested quantities or operators. Nevertheless, the ultimate theoretical values 
for such doings still remain to be seen, at least in nonperturbative regimes. 

On account of the above reasonings, it is definitely desirable or worthwhile to explore of any approach or 
procedure that would introduce no or as less as possible deformation of the bulk aspects of a quantum field 
theory and the associated infinities or divergences. For physicists, it is also more convenient to work with 
Feynman rules and Feynman diagrams that are physically intuitive and well defined at tree level in such a 
manner that the radiative corrections or loop diagrams could be computed without introducing too much 
mathematical auxiliaries. Actually, as will be demonstrated and argued below, there may exist an approach or 
strategy that seems to be free of the problems mentioned above or satisfy the requirements just enumerated. 
It is actually based on an elaboration of the effective field theory philosophy [2l| - |25| . introducing no artificial 
deformations, and more pleasantly, no ultraviolet infinities and the associated subtraction of infinities. Only 
local ambiguities for each loop integral may occur as a parametrization of our ignorance, which, as in the 
last step of renormalization, are to be fixed through appropriate boundary conditions. We feel all these 
virtues should make our strategy a simple and natural approach to start with. It has already been partially 
employed in a number of field theoretical issues[26-30]. 

From this report on, we wish to launch serial works to demonstrate, explain and further develop our 
strategy, with our ultimate goals targeting at a consistent and efficient program. In the course of these 
works, we will illustrate through concrete examples that many technical subtleties such as external momenta 
routing (or shift of integration variables) and overlapping divergent diagrams do not cause any trouble. From 
our demonstrations, one could find that our strategy could be well applied beyond the standard perturbative 
framework in terms of Feynman diagrams, see the applications in nonperturbative problems [29l |30| . unlike 
those principally devised for and limited to relativistic perturbative formulation. Another virtue of our 
approach is that we may examine important issues such as symmetry status of the radiative corrections 
(perturbative or nonperturbative) without being stuck to some specific prescription that might lead us to 
wrong conclusions or judgement [26- 28] . For perturbative amplitudes, our strategy could reproduce the 
EG approach's results. In fact, our physically motivated strategy share the following consensus with the 
mathematically motivated EG approach: the perturbative quantum field theories are not unambiguously 
defined in the short distances, some elementary Feynman or loop diagrams contain local ambiguous terms 
that could only be fixed through appropriate boundary conditions. In this report, we employ the simple 
and well-known objects in quantum electrodynamics (QED) for a clear and pedagogical illustration of the 
"naturalness" and simplicity of the principles and perspective that are adopted in our approach. In a sense, 
we wish to motivate an elaboration of the effective field theories philosophy with respect to renormalization 
of all quantum processes. 

The organization of this report is as follows: We describe our strategy in Sec. II along with some technical 
issues. In Sec. Ill, we compute the elementary one-loop amplitudes or Feynman graphs (that are ill defined in 
QED) using natural differential equations following from the existence of a complete underlying theory. The 
Ward identities between the vertex and self-energy diagrams will be demonstrated as a valid constraint and 



The reason is obvious: Were it found, it must at least be equivalent to the true underlying theory that is well-defined in every 
aspect, including the correct formulation of quantum gravity. 
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the four-photon (light by light scattering) diagrams will be shown to be finite as an additional consequence 
of gauge invariance. In Sec. IV, we consider how to impose appropriate boundary conditions to fix the 
ambiguities, how this step is linked to the conventional procedures, and related issues. We will also answer 
the origin of renormalization group equations in Sec. IV. The whole presentation is summarized in Sec. V, 
where some conceptual issues will also be addressed. 

II. DESCRIPTIONS OF OUR STRATEGY 

Our strategy is in fact to determine the ill-defined (or undefined) loop amplitudes through solving natural 
differential equations together with imposition of physical boundary conditions. In terms of Feynman dia- 
grams, this is to "calculate" divergent diagrams using convergent ones. To proceed, we denote a superficially 
divergent or ill-defined diagram and its underlying theory version with the same symbol T which should be 
a function of external momenta [p], couplings and masses [g] and the underlying parameters {a} (which will 
be usually hidden to avoid heavy formulae) that render the diagram finite or mathematically meaningful. 
As the underlying parameters {a} must be very small in "size", an ill-defined diagram should be determined 
through the "decoupling" limit of the version with the underlying parameters. Not knowing the details of 
the underlying theory, a natural and yet legitimate operation about the underlying theory version of an 
ill-defined diagram is to differentiate it with respect to external momenta [p] for enough times so that the 
resulting loop integration could be computed in terms of the conventional Feynman rules of quantum field 
theories 2 : 

^ +1 T(\p],{9})=f^)(\p],[9}), (1) 

with uir being the superficial divergence degree of diagram T and r("r+i) denoting the well-defined diagrams 
generated by the operation dp r+1 . Then we solve the above differential equation to arrive at a finite 
expression of T that contains a polynomial in terms of external momenta [p] with ambiguous coefficients. 
So, the general solution for the ill-defined diagram T reads: 

r([p], [g]) = f([p], [g]) + T^Qp], [g], [c]). (2) 

Here, [c] denotes the integration constants collectively. As last step of our strategy, the ambiguities are 
removed through imposing reasonable symmetries AND appropriate boundary conditions[31]- It is obvious 
that in our strategy, no artificial deformation is introduced at all. 

As our method rely on operations with respect to external momenta, there may be concerns about the 
issue of external momenta routing. Below we will show that: 

Theorem of Routing II. 1 The routing of external momenta in a loop diagram does not matter at all in 
our strategy, i.e., different routings will at most yield a reparametrization of the ambiguous polynomials. 

The routing of external momenta has been first addressed in Refs.JHJ and more than a decade ago. 
Here, we give a more concrete formulation of the remarks given there. There may be two kind of routings: 
(a) one keeps the distribution of external momenta at vertices intact; (b) the other amounts to relabeling 
the external momenta at vertices: [p] — > [p 1 ], or a redefinition of the external momenta at vertices. First, we 
prove it for one-loop amplitudes. 

Proof. Case (a): As the underlying theory version is well-defined, one could perform a integration 
variable transformation so that a different routing is exactly resulted, thus, 

r([p],[ 5 ];{a}) rtl =r([p],[ 5 ];W) rt2 , (3) 

with the underlying parameters {a} being now explicitly included for illustration. That is, this variable 
transformation should not alter the evaluation. Now apply the differentiation with respect to the external 
momenta &^ on Eq.(|3]), we end up with the following equality: 

02T([p],[g]; W) rtl =d^T([p],[g];{a}) rt2 . (4) 



2 In underlying theory perspective, the loop integration and "decoupling" limit (or low-energy projection) commute with each 
other on such differentiated diagrams, or that the "decoupling" limit can be performed before loop integration is done [2l!423| ]. 
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Taking the " low-energy" limit so that {a} become " decoupled" , we have: 

Limit w [T(\p], [g]; {a}) nl - T(\p], [g]; {<r}) rta ]} 

= d»][T(\p],[ g ]U-n\p},[9}U2}=0. (5) 

That means different routings differ at most by an polynomial r Poly that would be annihilated by 9£T : 

r p ^([p], [g]) ee T([p], [g]U - T(\p], [g]) rt2 , (6) 

which could well be absorbed into the ambiguous polynomials in T (to be fixed later through physical 
boundaries) as a reparametrization of the expression. 

Case (b): In this case the operation about external momenta should also not alter the amplitude provided 
the integration is done in underlying theory: 

r(b],[ 5 ];W) = r(b'],[ ff ];{a}). (7) 

Now apply the differentiation with respect to the external momenta [p]a that are in a subset shared by [p] 
and [p'} to Eq.0, we have, 

^jrM [<?]; M) - r([p'], [g]; M)] = 0. (8) 

Then we can take the "decoupling" limit and perform the loop integrations in the two routing, just as in 
case (a). Again, the net difference between the two routings is at most a polynomial with respect to [p]o that 
could again be absorbed into the ambiguous part, after we return to the original label of external momenta, 
the same expression for T should be restored. 

The multi-loop diagrams could be treated in similar fashion 2 lj-|25j: First, differentiate the diagram with 
respect to the momenta that are external to the overall diagram to annihilate the overall divergences; Second, 
consider each individual sub-diagram that is divergent and differentiate it with respect to the momenta that 
are external to this sub-diagram to annihilate the "overall" divergences of the sub-diagram; Third, continue 
this operation till the smallest sub-diagrams are thus differentiated; Fourth, carry out all the resulting loop 
integrations and integrate back indefinitely with respect to the corresponding momenta that are external 
to the corresponding loops; The final outcome will be a finite expressions in terms of momenta and masses 
external to the whole diagram, containing ambiguities associated each individual divergent or ill-defined 
loops. This operation naturally dissolves the overlapping divergences that are notoriously difficult to deal 
with in conventional treatments, thanks to the work by Caswell and Kennedy (32j. Since the treatment should 
be done loop by loop, thus, the theorem of routing applies to each loop integration as well and finally to the 
whole diagram. Q.E.D. 

This theorem would allow us to choose any routing we like or that is convenient. Obviously, simple routing 
of external momenta would yield simplicity in calculations. 



III. QED AT ONE-LOOP LEVEL 



We will work with the following standard covariant gauge QED in 3+1-dimensional spacetime: 

£qed = # - m)V - \F^F^ - ^(d,A») 2 , (9) 

with £ being the dimensionless gauge parameter. Our metric convention reads 

9^ = 9^ = (+,---)• (10) 

At one-loop level, only the following five elementary vertex functions are superficially UV divergent and 
hence have to be renormalized: the self-energy E, the electron-photon vertex A M , the vacuum polarization 
ILj„, the three- and four-point photon vertices. The Feynman diagrams for the first three vertices are given 
in Fig.l. The three-point photon vertex is zero due to Furry 's theorem, while, as will be shown below, the 
four point photon vertex is actually definite in our approach due to gauge invariance. 
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We first consider the vertex functions listed in Fig.l whose Feynman integrals read, 
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where p, p and q denote the external momenta for electron and photon respectively. The superficial divergence 
degrees for the three diagrams are well-known: u;(£W) = 1, lu(A^) = 0, w(H$) = 2. Note the routing 
of the external momenta are so chosen that they flow through the fermionic internal lines, which will yield 
convenience for calculations and a clear diagrammatic interpretation of the operation associated with the 
differentiation with respect to external momenta, see below. We will consider other routings for the self- 
energy and the vertex diagrams in the Appendix A. 

A. Differential equations for divergent integrals 

Below, we treat these superficially divergent integrals with the strategy described above, according to 
which the three differential equations are, 



d Pa A^(p,p) = A^ a ( P ,p), 

d qa d qfi d q n^(q)=U^(q), 



(14) 
(15) 
(16) 



where the right hand sides (RHS) of these differential equations are convergent Feynman diagrams, their 
corresponding loop integrals are listed in Appendix A. Actually, in terms of Feynman diagram the RHS of 
Eqs. (|14ll5ll6[) are just obtained through insertion of photon probes with zero momentum, provided that the 
external momentum flows through internal fermion lines, cf. Fig. 2 for an example. 

The loop integration on the RHS of Eqs. (|14ll5ll6[) can now be straightforwardly carried out: 
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Here symmetrization about Greek indices grouped in the bracket {•••}) is frequently used for brevity. The 
three definite functions then determine the original loop amplitudes up to certain degree of ambiguity. Next, 
we proceed to solving the differential equations using the above definite functions. 



B. Solutions 



Evidently, the direct integrating back with respect to external momenta would be cumbersome. Instead of 
doing so, we may first extract the differentiation with respect to external momenta in the following manner, 
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with [C...] (C M „ is traceless, = 0) being the corresponding integration constants that are ambiguous at 
present stage. They are uniquely defined only in the underlying formulation that is well defined at short 
distances, but in the effective theory, i.e., QED, they are free constants to be determined through appropriate 
boundary conditions [31 1. 

Actually, for photon vacuum polarization, Lorentz invariance could remove the constants C^C^. The 
presence of Cn ; o would violate gauge invariance, i.e, making the photon massive. To ensure gauge invariance, 
we must impose C a + C 'a — Cn ; o — 0. For self-energy, both and C m preserve Lorentz invariance, but are 
gauge dependent. In our previous reports, the constants [C...] are termed as "agent" constants as they should 
be consequences of taking the "decoupling" limit of {cr}. To further determine the residual constants, we 
must resort to more constraints from symmetries and finally from physical boundary conditions, which will 
be pursued in the following section. Before imposing boundary conditions, the net outcomes of conventional 
renormalization programs are also ambiguities, the true reflection of ill-dcfinedness. The same is also true in 
EG approach. However, our strategy leads to the same outcome with neither complicated (also unnecessary) 
manipulation of infinities nor modification of Feynman rules at loop levels, hence, it is simple. So, the real 
problem is ambiguity (rather than divergence) to be fixed through appropriate boundary conditions. 

Here, some more technical remarks are in order: (1) In the differential equation approach, any spurious 
violation of symmetries due to local ambiguities in the loop amplitudes could be easily removed. Therefore, 
the physical breakdown of symmetries (i.e., anomalies) must come from definite properties of the loop 
diagrams. Actually, it is shown in our previous studies that chiral and trace anomalies [33Tj4]]| d o originate 
from the existence of definite or nonlocal terms like (k^k„ /k 2 )T in certain divergent integrals 3 [26l I27I. |42|. |43| . 
After contraction with /c M or g^ v , they give rise to local operators: (a) fc^ • (k^k v /k 2 )T = k v T (chiral 
anomaly) 26, 27, 42], (b) g^ u ■ (fc M fc„/fc 2 )f — f (trace anomaly) (421 |43| . In contrast to spontaneous breaking 
of global symmetries where massless poles appear at tree level as consequences, certain massless poles show 
up in loop amplitudes as origins of the quantum mechanical violation of canonical symmetries. (2) Amusingly, 
this perspective of anomalies also imply that there may exist properties in divergent integrals that are both 
definite (nonlocal) and unexpected from canonical deductions. Therefore, such divergent diagrams contain 
more about quantum field theories [HI, EH- These properties might be viewed as indirect but definite links to 
underlying theory. (3) The differential equations used here also imply some structural relations between the 
Feynman diagrams. Actually, in QED, as there is no pure gauge boson loops, all the differentiations with 
respect to external momenta effectively amount to inserting the elementary QED vertex, or external photon 
with zero momentum. The final result amounts to the Ward identity in differential form. It will be interesting 
to see what will happen in more complicated theories with pure gauge boson loops, like non-Abelian theories. 
(4) Technically, one could not obtain a polynomial/local term from convolution of non-polynomial/nonlocal 
functions. A mathematically natural way to yield a polynomial/local term is to perform certain limit 
operation. To our interest, this mathematical scenario corroborates the underlying theory scenario: The 
computation should be followed by the " low-energy" limit due to the wide separation of scales, then some local 
terms arise from the "low-energy" limit, indicating the existence of underlying structures. Conventionally, 
a regularization is introduced as a rough substitute of the underlying structures and is taken to zero in the 
final stage after removing the "dusts" (infinities) brought about by the regularization. In our strategy, we 
simply appreciate the existence of the underlying structures and then solve the differential equations that 
must be satisfied by the limiting objects. 



C. Comparing with Dimensional Regularization 



Now we list out the results from dimensional scheme for a comparison. Our conclusion is that at one-loop 
level the only difference between the two approaches lies in the local part, which is in fact prescription depen- 
dent. The definite (nonlocal) part must be prescription independent and hence physical. The two approaches 
are equivalent if one takes the local part in dimensional scheme as ambiguous. From this equivalence one 
could devise a simple method for computation: do the calculations in dimensional regularization, then re- 
place the local part with ambiguous polynomials of the corresponding degree. However, this equivalence is 



3 Here, k denotes the external momentum at the corresponding vertex, while f the additional polynomial factors in terms of 
external momenta at other vertices. 
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only valid after sub-divergences are subtracted in case of multi-loop diagrams [46 
The results from dimensional scheme read, 
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We note that in terms of the above parametrization, the following correspondence between the agent constants 
and the local constants in dimensional regularization could be easily established: 
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Obviously, the dimensional regularized results exactly agree with ours obtained from Eqs. (|14|15|16p after 
replacing the divergent constants by the agent constants following the correspondence (|30[) . In other words, 
the dimensional regularization results (after or before subtraction) could be seen as a particular solution to 
our differential equations at one-loop level. At multi-loop level, the dimensional results will also satisfy such 
differential equations loop by loop after all the corresponding sub-divergences are removed (46|. 



D. Reducing ambiguities with Ward identities 



Since gauge invariance is encoded in the Ward identities among various vertices, it is natural to ask if these 
identities are satisfied for the loop amplitudes computed in our approach. Specifically, we wish to examine 
the status of the following Ward identity for self-energy and vertex function, 
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To proceed, we make use of the results in dimensional regularization and look at the Feynman gauge 
components, i.e., those remain when £ = 1, the rest components will be delegated to Appendix B. Carrying 
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(5^ - 2m) + i> 



6 

2 



4m 



(32) 



1) In ■ 



5-1 



(47T) 



2 5 In 



5 y j p 

Cj_(e)||5=i + 3 - In m 2 + 5 + (5 2 - 1) In ; 

5-1 . A , „ , _ , ] 

P 2 



7m 



7p 



1(5/*- 2m) + i> 



(33) 



with 5 = 2pj-. The constants C^(e), C±(e) have been defined in the correspondence (|30l) . Now differentiating 
EW(p; e) with respect to p M , it is trivial to verify the Ward identity (|3ip in dimensional regularization: 

5-1" 



(47T) 



2 51n 



7-B 

- 1 



1 - In 



1 



47T^ 2 

(5/* — 2m 



+ 5 + (5 2 - l)ln 



7^ 



Pm 

„2 



(34) 



= iAW(p,p; e). 

It is easy to see that in dimensional regularization this Ward identity is ensured by the following relation, 

In general gauge, we have (Cf. Appendix B): 

Ci ,(e) + C ± (e) = ^lt. (35) 



Evidently, replacing C^(e), Cj_(e) with C^,C±, one enters our approach. Thus the Ward identity (pTTj) 
could hold in our approach only when the agent constants and C± satisfy the same relation as that 
between C^(e) and C±(e): 



Cj, + Cj 



11^ 



(36) 



That means, gauge invariance constrains the agent constants C'^p and C± according to Eg. (1361) . Obviously, 
such constraints are just welcome in our approach to reduce the ambiguities. Of course, at higher orders, 
such constraints must be also consistently imposed. At one-loop level, there is in principle no obstacle to 
remove the local violations of any symmetries. As removal of the local ambiguities never affects the nonlocal 
and hence definite part, any physical anomaly must be originated from definite and hence prescription- 
independent sources. As mentioned above, at least for chiral and trace anomalies, this is indeed the case: 
a type of mass-independent and definite terms is the true source of quantum mechanical violations of the 
canonical chiral and scale transformations (2(| [2?], HU, [43|. 

Therefore, due to Lorentz and gauge invariance, we are left with three ambiguous constants to fix: C^,, 
C m and Ca- This statement is, however, a little bit cursory before the finiteness of the rest superficially 
divergent diagrams, i.e., the three- and four-photon vertex functions, is established. Below we turn to 
these diagrams. Since the three-photon vertex does not contribute to cross sections or physical observables 
according to Furry's theorem, we only need to verify the finiteness of the four-photon vertex, that is, there 
is no ultraviolet ambiguity in this vertex function which is again ensured by gauge invariance. 
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E. Finiteness of the four-photon vertex 



To verify that the four-photon vertex is free of ambiguity, it suffices to show that the most divergent piece 
in each individual diagram cancel out against each other due to gauge invariance. 

The external momenta and polarizations of the four external photon are arranged in the following manner: 
(Pi> A 4 )) (P2, v), (p3, p), (p4, a) with p4 = — pi — P2 — P3 due to conservation of the four dimensional momentum, 
i.e., pi always enters at vertex /i, p2 at v, etc. The four-photon vertex consists of six inequivalent individual 
Feynman diagrams that differ from each other by the arrangement of the relative positions of the external 
photons along the internal fermion loop. It is convenient to fix the position of the first photon of (pi, /z) and 
let the other photons move to obtain the rest inequivalent loop diagrams. The notations for the four-photon 
vertex and the six individual diagrams are as follows, 

f pi>pa\Pl: P2: P3: P4} — ^ l;pi> pa {pi,P2,P3,Pi) + r 2 (Pl,P3,P2,P4) 

+?3- tJ .vap(Pl,P2,PA,Pz) + ^4;papiy(Pl,P4,P3,P2) 

+^h:,paup{Pl,P4,P2,P3) + T 6]tiptrv (jpt , p 3 , p 4 , P2 ) ■ (37) 

The Feynman integral for Y\.^ vpa reads, 

4 f —d 4 l ( i i 
J (27T) 4 \/, — m + ie /l+pi—m + ie 

X H+i>i+j)2 - m + ie lp /I TJi +i>2 +h - m + i€ la ) ' ^ 

By rearranging the external momenta and the elementary vertices, one could readily obtain the integrals for 
the rest five diagrams. Since the superficial divergence degree of such diagrams is 0, then it suffices to verify 
that the logarithmic divergence is absent in the sum of the six diagrams in order that gauge invariance is 
preserved. Obviously, it suffices to do computation with Y\. pvpa as others could be obtained by permutations 
of the vertices (pi, //), (p2, v\ ■ ■ ■. 

Differentiating the integral in Eg . ([551) with respect to pi- a , the resulting integral becomes definite and 
gives, 

8ie 4 f 1 f x f v 

dpi-JTi-uvpa = ^ T\;u,vpo dx dy dzd Pl a In A x + d Pl .JT Ufl „ p<7 , (39) 
(4tt)^ Jo Jo Jo ' 

Tl- : p,i/ pC r = Qpvdpo ^9pp9va QpaQvpi (40) 

Ai = m 2 - (x - x 2 )pl - (y - y 2 )p\ - (z - z 2 )p\ 

-2y(l - x)pip2 - 2z(l - x)pip 3 - 2z(l - y)p2P3, (41) 

where the Y\-^ vpa denotes the definite parts after integrating back with respect to pi ;a and is not concerned 
here. Next, integrating back with respect to pi ;Q , the diagram T\-^ vpa could be parametrized as follows: 

8ie 4 f 1 f x f v Ai 

iTi-uupa = ,, Tz Ti-uvpa I dx I dy I dzln — \-T\-u Vp<T , (42) 

(4tt) z J J q J q Gi;4 7 

where an integration constant Ci ; 4 7 is shown explicitly, which is exactly the ambiguity that corresponds to 
the logarithmic divergence of the integral in Eg. (1381) . The results for other diagrams are given in Appendix 
C. For later convenience, we integrate by parts with respect to the Feynman parametric integration to arrive 
at the following expressions for Y\. illvp(T , 

l-iiv per — n i a \0 -^"ljAlfPCT kl ^Muvpcr: (43) 

where the definite terms arising from this operation have been absorbed into Ti- pup<T and collectively denoted 
as ri;p„ pcr , cf. Appendix D. 

Now, in order that the ill-defined piece in the full amplitude vanish 



T II - J!f_ 

«■ M"pfllill-def ined — 3(47r) 



Jfe=l 



0, (44) 
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with Ylt=i Tk;nvpa = 0, the constant Cfc ; 4 7 must be identical for each k (say, Ck-A-y — C < 4 7 ,Vfc). Otherwise, 
gauge invariance will be violated. Thus, similar as the case of vacuum polarization for two-photon vertex, 
gauge invariance also remove the ambiguities in the four-photon vertex function. 

In fact, after setting all the external momenta to zero, we will end up with the same ill-defined integral in 
all the six diagrams: 

dH IJblJd i r m 2 

(2tt)* (Z*- m 2)4 = 2A(^ [9ab9cd + "- ]ln cT 1 - (45) 

Thus it is natural to deem that the constants Cfc ; 4 7 are identical in any reasonable prescription. 

In dimensional regularization, one would find the following correspondence in each of the six diagrams: 

m 2 i , m 2 

In — -o- -e + 7b + In ■ 



Cfc;4 7 ' 47r^ 2 ' 

which means the full four-photon vertex is finite in dimensional regularization due to the mechanism demon- 
strated above. 



F. Brief summary of differential equation approach to QED at one-loop level 



Now, we have treated all the superficially divergent or ill-defined loop diagrams of QED at one-loop level, 
without resorting to any artificial deformation of spacetime and symmetries or any form of cutoffs. At one- 
loop level, only three local ambiguities remain to be fixed in QED: C$, C m and Ca that appear in self-energy 
S, vertex and vacuum polarization 11^, respectively. This is a natural conclusion following from Lorentz 
and gauge invariance and the existence of a complete underlying theory. 

In particular, from our treatment, one could naturally infer that, any fundamental theory that is well de- 
fined and underlies the presently known QED must also yield the same or equivalent functional expressions as 
given here after taking the "decoupling" limit, the only difference is that the unknown agent constants could 
be unambiguously computed in such underlying theory. Therefore, in the underlying theory perspective, the 
QED processes defined by the elementary Feynman diagrams (in covariant gauge) represented by Fig.l must 
take the following form, 



A^(P,P) 



in#(g) 



ie 



(4*0- uo 



dz 



[(l-2z-^+(3 + eHln m2 |2 ^ 



(1 - 04P 2 - m 2 )i> 



m 2 — zp 2 



(47T) 



1 rl-y 

dy I dz 



Mud 

+ (%i> + c° n m j , 

2 7M lnA + ^ + (l-0(l-y-z) 



i 



4 + ne o , 

—7 — <Vpf 



8ie (9/j.vq - g^v) 
(4tt) 2 



, . _ TTh z(l - z)q 2 n ' 
dzz(l - z) In \ '— + c° A 



o 



Mud 



(46) 



(47) 



(48) 



where we have explicitly "separated out" a dimensional constant ^ud from all the constants C... to balance 
the dimension in the logarithmic parts, the residual parts are purely dimensionlcss numbers, each denoted 
as c ... This way of parametrizing the constants will be more precise and useful for the following discussions. 

We should also note that, in the underlying theory perspective, all the constants or parameters that 
appear in the Lagrangian should be defined from the "decoupling" limit of the underlying theory also as 
" agents" . In our view, these " agents" in lagrangian are only elementary from the perspective of (effective) 
field theories, not necessary so in the perspective of underlying theory. As the agent constants appear in the 
local part of the elementary vertex functions within perturbative framework, the Lagrangian constants would 
mix with these local terms containing the agent constants in the corresponding vertex functions. This is in 
fact the origin of renormalization group equations in our approach, without resorting to the conventional 
renormalization theory, see, Sec. IV for a brief discussion. 
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IV. REMOVING AMBIGUITIES AND " RENORMALIZ ATION" SCHEME 

Now we turn to determining the ambiguous agent constants. Obviously, the conventional renormalization 
schemes would be reproduced provided that the ambiguous constants as well as mass and couplings are fixed 
through the conventional renormalization conditions. An appropriate choice of the whole set {m, e; /x, c ..} 
defines a renormalization prescription or scheme. So, our simple strategy is compatible with the conventional 
renormalization programs. In principle, any scheme can be employed provided the same scattering matrix 
elements or the same physical observables could be obtained. This freedom just leads to the renormalization 
group in the general sense of Stuckelberg-Peterman[47|. Fixing the dimensionless constants in a scheme and 
letting the scale /irjD vary, one could obtain the renormalization group equations in the narrow sense. 

To proceed, the following parametrization would be used useful, 

u(p)K^u{p) 

n#(g) 

F^(k 2 ) 



G {1 \k 2 ) 

n<V) 

5 

Note that the vertex function is sandwiched between on-shell states [u(p),u(p)] of electron. Obviously, the 
scalar factors A^ x \ and IlW contain ultra-violet ambiguities (depending on c°,c^,c^ and c° A , 

respectively). Moreover, the first three factors are also gauge dependent and infrared singular when the 
external momenta are on mass shell. In we have explicitly put in a photon mass squared (A 7 ) to 

regulate the infrared singularity, with S denoting the terms proportional to 1 — £. In contrast, the scalar 
factor G^ 1 ) is both ultra-violet and infrared finite and also gauge independent. Its value at zero momentum 
transfer (k 2 = 0) gives the well-known anomalous magnetic moment of electron at one-loop level: 

G<..(0) . £. 

A. Various boundary conditions for vertex function 

Now we impose various conditions on the elementary vertices (here, T,^\ A^ 10 , nj^) 4 . In doing so, both 
the tree level parameters (m, e) and the agent constants must be understood as scheme-dependent so that 
the observables computed using these parameters and vertices are physical and independent of schemes. 
The lagrangian or tree level mass or coupling needs not be just the physical one, which may differ by a 



A {l) {p 2 )i> + B {1 \p 2 )m, 

u(phMP)F {1) (k 2 ) + ^%^) G a) (fc2)! 



in ' 

-n (1 V)Gv<r- <i„h, 

a ' 

ill 



2m 



a 

47T 

a 
-in 



£( {6 2 -l)ln S —± +8-ln^-) 
X 6 Mud / 

A — 1 777 ^ 

(3 + ( (5-l)ln— + l-ln — 
6 Mud 



2arccoth[6'fe] 



2 + fy 















■k) 




+ 2 













-S(fc 2 ;l-C)-ln 



Mud 



a 
3^ 

9 

m 



q {9 2 q - 3) arccoth[0 g ] - Q\ - In 



Mud 



?-6c° 



= —o, k = p-p, e x = \ i - 



9 ' f fl "X — \/ * 9 



4m 2 



(49) 
(50) 
(51) 
(52) 

(53) 



(54) 
(55) 
(56) 
(57) 



4 They could be found in many textbooks for field theory, see, e.g. Ref . |4S( 
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" (re)normalization" constant in the conventional terminology. In our strategy, such " (re)normalization" 
constants could also be introduced but are nevertheless ultra-violet finite. Actually, as is already pointed out 
in Ref.[49(: "the renormalization of mass and field has nothing directly to do with the presence of infinities, 
and would be necessary even in a theory in which all momentum space integrals are convergent" . 



1. On-shell conditions 



First, let us consider the so-called on-shell conditions in Feynman gauge, which read, 
£(% =m = 0, ^£(% =m = 0, AW||^ =m = 0, n«|| 92=0 = o. 
Using the expressions given above, we could find that: 



(58) 



a 
An 



— — -31n 



Mud 



a 
An 



= — In 



2 ^ 
A, 



2 

Mud 



-21n^---d 



id A 



a 
An 



In 



Mud 



a 
Yn 



= ^ -In 



Mud 



-6c° A 



(59) 



where 



= 5/2 has been used. Note that all the infrared singularities have been regularized with a 
photon mass, and m now denotes the physical mass for electron. According to conventional programs, the 
coupling e or a = e 2 /(An) appearing in these conditions is understood to be the physical one. 

Then the on-shell conditions could be fulfilled provided the constants [c ..] take the following values: 



In 



Mud 



2 ~ i A-y 

2 In — 1 

m 



7 



- 2 2 ' 



-4 In 



m 

Mud 



C A 



1 772 

6 Mud 



In this manner, all the ambiguities are gone in the vertices, including the scale /i\jd- 



a£Hp 2 ) 


a 

~ An 


b£Hp 2 ) 


a 

~ An 




a 

~ An 
-4- 


iW) 


a 

= 3^ 



(5 2 - 1) In 



5-1 



5 - 3 - 2 In 



A 



4 ( < S-l)l n __i-2-21n^ 



2arccoth[6'fc] 
A. 



In ^ 



21n^ 

772 z 



(9 2 - 3) arccoth[6> 9 ] 



(60) 

(61) 
(62) 

(63) 
(64) 



The on-shell scheme works for field theories without unstable fields as elementary quantum degrees of free- 
dom. The tree level mass is identified as the physical pole mass. But for c omp licated theories like electroweak 
with unstable sectors, the on-shell scheme runs into trouble, see Ref . |50h59| . 



2. "Minimal" conditions 



In analogy to the famous "Minimal Subtraction" , one may define a scheme where [c ..] are removed so that 
only the ~ In \j? pieces are left over, which will lead to the well-known mass-independent running behavior. 
We will denote this scheme as " Ms" , which may also be interpreted as " Minimal specification" . 
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Since A and can not be zero at the same time due to gauge invariance, we choose the following: 



r° 4- - 



C A 



0. 



(65) 



Then the prescription-dependent scalar factors would take the following neater forms (in Feynman gauge): 



Klip 2 ) 



O 2 ) 
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l)ln 



(5-1 



5 -In 



Mm s 
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7T 

a 
4tt 



/ r > . 5 1 „ , to 
{5 - 1) In — h 1 - In 



Mm s 



2arccoth[6'fc] 
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-1 - In 



n (D f 2) = a_ 
M ° [q ' ~ 3tT 



9 q (d\ - 3) arccoth[6» g ] - . 



h + 



8 , to 
3 Vlis 



(66) 
(67) 

(68) 
(69) 



Here the subscript " Ms" for mass m and coupling e are omitted to avoid heavy symbolism and the scale jUud 
is now replaced by ji Ms ■ Of course, the concrete values for to and e in this scheme need to be determined in 
terms of physical values of mass using appropriate physical observables 31] . In this prescription, the scalar 
factors A and B at general off-shell momenta are not beset with infrared singularities in contrast to the 
on-shell prescription. 



3. Boundary conditions at vanishing off-shell momenta 



One could also impose boundary conditions at off-shell momenta. Among these off-shell prescriptions, the 
zero momenta case is much simpler for QED with massive fermion. In Feynman gauge, we have: 
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Then the scalar factors read: 
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(70) 
(71) 

(72) 
(73) 
(74) 
(75) 
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The nice point of off-shell conditions at zero momenta (hence the subscript "zm") is that, the infrared singu- 
larities do not show up in the scalar factors at general off-shell momenta, just like in the "Ms" prescription. 
Such conditions possess an intuitive interpretation: at zero momenta, quantum fluctuations are suppressed, 
and an electron would serve as a classical static source of electromagnetic field, with the tree level coupling 
" e" being pertinent to static or Coulomb charge somehow. 

4- Physical boundary conditions in terms of observables 

In principle, boundary conditions may be imposed upon many other field-theoretical objects that are 
constructed with elementary vertex functions or Green function. It would be better to work with the 
observables that could be readily confronted with physical data and depend on the elementary parameters 
(masses, couplings, and agent constants) in an as simple as possible manner, sparing the intermediate 
renormalization. 

For example, one may proceed as below in QED: First, physical values are assigned to the Lagrangian 
parameters (masses and couplings) , the agent constants are treated as unknown or arbitrary; Second, appro- 
priate observables are selected that could be readily expressed in terms of Feynman diagrams that definitely 
contain the agent constants; Finally, the agent constants could be fixed in terms of the physical parameters 
and data of these observables. In formulae, 

(1) {m,e; [C...]} => {m (phy) , e (phy) ; [C (phy) ]}, 

(2) Q n (m(P hy ), e(P hy ); [C. (phy) ]) = Q< data) , n = 1,2, ■ • • 

=>• C (phy) = /... ([Qi data) ]; m (phy) , e (phy) ) . (76) 

The set of observables to be employed for this purpose should fulfill some natural requirements such as 
simplicity and analyticity with respect to their dependence upon mass, coupling and the agent constants, well- 
definedness in the infrared, etc. In principle, even theories like quantum chromodynamics (QCD) that are 
elusive of physical pole masses of elementary quantum fields can be treated in this manner, where the physical 
contents of the lagrangian or tree level masses and couplings becomes a nontrivial issue. Such imposition 
of physical boundary conditions is general and natural from underlying theory perspective. Actually, this is 
exactly what is done in literature for the determination of the "physical" coupling of strong interaction, see, 
e.g., Ref.[ll|, chapter 12. 

5. A remark on boundary conditions 

In order that the same observables are obtained, different choices of "intermediate" boundary conditions 
or prescriptions for vertex functions must be related to each other somehow, or more specifically, the trans- 
formations across different prescriptions defined by the boundary conditions should in principle be feasible, 
at least within perturbation theory. In conventional programs, this is the so-called scheme dependence in 
perturbation theory, where different renormalization conditions or prescriptions, differ by a set of finite 
renormalization constants, and the observables computed in different prescriptions should agree with each 
other to the perturbation order computed. This is in fact encoded in the renormalization group equation 
in the sense of Stiickelberg-Peterman as mentioned in the very beginning of this section. Furthermore, a 
general prescription would also contain an arbitrary scale as " renormalization point" . The variation of this 
arbitrary scale within the specified prescription should not affect observables, which would in effect lead to 
renormalization group equations that are themselves prescription- or scheme-dependent. 

The above statements make perfect sense within perturbative framework. However, things might become 
complicated in nonperturbative contexts. As an simple example, we refer to our previous work[28|, where 
it is shown that some scheme in incompatible with dynamical symmetry breaking of massless A</> 4 theory 
in nonperturbative context. Strictly speaking, any intermediate boundary conditions or prescriptions to be 
employed must be in the same "orbit" as the physical boundaries. This in turn requires a complete solution 
of the full theory of quantum fields, which is presently out of our reach. So, scheme dependence is a real 
issue for quantum field theories beyond perturbative regime. 
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6. Influences of infrared divergences 

It is well-known that QED is beset with infrared divergences. For on-shell scheme, infrared divergences 
show up already in the boundary conditions for the loop amplitudes, as is evident above. Simply speaking, 
infrared singularity arises because Fock states could not work well in presence of massless particles. Thus 
introducing certain extent of coherence, the infrared divergences in QED could be removed by use of Bloch- 
Nordsieck theorem [60|. which has been described in detail in many standard field theory textbooks, e.g., 
[3ll |48| . l49j . At one loop level, one regularize the infrared singularities in the loop diagrams somehow (here 
in this report, in terms of a fictitious photon mass) and they will cancel out against those from real diagrams 
in physical observables. 

Of course, for more complicated non-Abelian gauge theories, severe infrared divergences are present and 
their treatments are more involved. Nevertheless, as our approach do not alter anything in the infrared, so all 
the conventional methods for dealing with infrared singularities or similar singularities except the ultraviolet 
ones may well apply. As is clear from our above discussions, one may work with a set of renormalization 
conditions that are not plagued by infrared singularities to avoid further complications, this could be done 
with conditions or "boundaries" defined at off-shell momenta where everything is infrared finite (6ll. |62|. e.g., 
at zero four-momenta as illustrated above. 



B. Lagrangian perspective 

In this subsection, we reexamine the issue in terms of lagrangian or action, which allows us to examine 
a number of issues in a more transparent manner. Such effort is also helpful to resolve some conceptual 
difficulties that have been "afflicting" the conventional practices. 

1. Brief review of conventional programs 
The conventional program of renormalization starts with a bare lagrangian: 

£(i> B ,A B ;m B ,e B ) = $ B (ipB - m B )vb B - -F B ..^F^ - -L(5^) 2 , (77) 

4 Zl; B 

with D B = — ie B A B . Then all the 1PI vertices (r B ([p l7 ■ • • ,p n ], m B , e B )) are renormalized through the 
following renormalization rescaling, 

ipB = Z x >\, A B = Z]( 2 A^ m B = f^m, e B = (Z A )~ l / 2 e, ( B = Z A £, (78) 
T ' Z^ 

and 

£(<ip B ,A B ; m B ,e B ) = Z^tf)^ - Z m m^ - ^-F^F^ - ^(d^) 2 

= C(i/>, A";m, e) + A£(^, A"; Z^, - 1, Z m -1,Z A - 1), (79) 

so that 

T B l) ([pi, ■ ■ ■ ,p n ],m B ,e B ) = Z^ /2 Z A nA/2 T^([ Pl , • • • , Pn ), m, e). (80) 

Here, AC serves to provide all the necessary counter-terms for canceling out any divergences in the loop 
diagrams. As the BPHZ algorithm is equivalent to the above program, below, we will refer to both as 
conventional approaches or programs. One caveat is that the formalism for such programs is only established 
within the realm of perturbation and hence must be so implemented. 

2. Underlying theory perspective 

Now, let us reexamine the same issue from underlying theory perspective. It is convenient to work with 
the path integral formulation. In terms of a underlying theory description, the generating functional for the 
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QED processes would formally take the following form: 

Z (J»,V,V;W}) = Jd{i[<j> VT ]exp^i Jd i x [Ajt([0ut], M) 

+JMM) + 0^{{a})n + vO^({<j})} } , (81) 

where [</>ut] is used to label all the necessary degrees of freedom, and 0^,0^,0^ denotes the photon and 
electron field parameters in the UT description which may well be composite ones. 

At "lower" scales where QED become prominent or effective processes, all the typical higher energy modes 
are actually integrated out, resulting in a well-defined functional in terms of dominant or effective degrees, 

Z(J fi , V ,fj;{a}) = JdA^d^exp^i Jd A x [C cS (^, A^; {a}) 

+fa + fri>]}, (82) 

with J d 4 xC e s(ip, A^; {cr}) being the effective action generated from the integrating-out. Taking the "decou- 
pling limit" (Vle) naively on the effective lagrangian would yield to a classical lagrangian of QED, 

A E [£eff(V>, A^- {a})} = £ Q ed(<A, AJ, (83) 

which in turn leads to ill-defined path integral as the limit operation and the functional integration does not 
commute. Obviously, the details taken away with the limit operation, 

5C(iP, A^; {cr}) = £ eff (V, A„; {cr}) - V LE [C cS (^, A^ {a})} 

= £ eff (V',A /J ;{<7})-£QED(V',AJ' 

is just what we are missing for a well defined path integral. Therefore, we must include it into the path 
integral, 

Z(J^,r],fj;{a})= JdA^d^exp^i J d 4 z [£q ED (V>, \) + SC(i>, A^; {cr}) 

+j„a» + $ v + ffti>]}. (84) 

Evidently, in conventional programs, SC is implemented through regularization and subtraction that are 
encoded in AC. For certain interactions, the job of AC could be "autonomously" done using operators 
appearing in the classical lagrangian, leading in effect to the "renormalization" of these operators. Such 
autonomous cases are conventionally termed as renormalizable, while the rest as un-renormalizable. As a 
"realization" or substitute of SC, AC must also take the responsibility to clear all the side effects like ultra- 
violet infinities associated with a regularization. Only in this sense, the manipulation of infinities (through 
counter terms) may be justifiable, and the bare lagrangian may serve as a symbolic and rough representation 
of the underlying theory: C c s{ip, A^; {cr}), which definitely contains the sophisticated underlying quantum 
details that could not be simply described in terms of "bare" objects that are originally coined in classical 
field theory. 

As we could at best obtain the parametric form of the "decoupling" effects of {cr}, there may be certain 
arbitrariness in our choice of £qed = PleICcs] (that is actually unknown to us yet) as our "starting" 
lagrangian for calculations or the in the separation of it out of C e s'- 

£ e ff = <£qed + SC(- ■ ■ ; {cr}) = £qed + 6£(- ■ ■ ; {cr}), 

where £qed and £qed would at most differ by a series of local operators (for renormalizable theories, by a 
finite "renormalization" of the lagrangian) that could be absorbed into SC(- ■ ■ ; {cr}): 

5C(- ■ ■ ; {cr}) = <5 var £QED + 8C(- ■■;{<?}), ^var^QED = £qED - ^QED- 

Obviously, any sensible prescription of V-us\C e s\ should necessarily include a reference scale that are widely 
separated from the underlying ones. This natural appearance of a reference scale just corresponds to and 
hence remove the mystery around the dimensional transmutation phenomenon in the conventional programs. 
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C. Origin of renormalization group equations 



We note that renormalization group (RG) and RG equations (RGE) appear in terms of intermediate 
renormalization reparametrization, for which a lucid discussion can be found in Ref. (o3| . If we could proceed 
with all parameters being physical, there seems to be no room for it. However, as discussed above, even 
starting with a physical parametrization, one could still arrive at an arbitrary prescription, provided the 
prescription is related to the physical one via proper transformations. In this short subsection, we will not 
repeat the derivation of RGE in a specific scheme defined by the corresponding boundary conditions as is 
done in conventional programs, instead, we briefly show that the RGE like equations could be generically 
derived as a "decoupling" theorem in terms of underlying structures even if the parameters are determined 
by physical boundary conditions [63 - [66| . 

First let us look at the scaling law of QED in terms of underlying theory. For any 1PI n-point vertex 
function T^, the scaling law would read, 

r (fl) ( [Ap] , [Am, e] ; { \ d ° a}) = A d r<~> r W ( [p] , [m, e] ; {a}) , (85) 

with d... denoting the canonical scaling dimension (equals to the mass dimension) of the corresponding 
parameter. The differential form of this scaling law reads, 

1 1> ' dp + m dm + J2 d ° ad ° ~ d r(-> I r (n) (M, [m, e]; {a}) = 0. (86) 

K P a ) 

Now applying the "decoupling" limit to this differential equation, we have, 



{& ' d P + md ™ + Y, d ^ d c - dr W | rW([p], [m, e]; [Q]) = 0, 



(87) 



with 



denoting the contributions arising from the "decoupling" limit operation. 

Note that, in whatever prescription for the agent constants [Cj], the scaling of these constants would just 
elicit insertion of appropriate local operators with respect to each loop where agent constants show up: 

Y,d Ci adcT^ =J2 S oJoT {n \ (89) 

i a 

where Io a denotes the insertion of a local operator O a that corresponds to a vertex whose 1PI loop corrections 
contain local ambiguities with So a being the associated coefficients. At least at one-loop level in QED, only 
logarithmic ambiguities will be really contributing to these " decoupling" effects that are anomalous in terms 
of the canonical contents of quantum field theories. Thus, So a is the primitive "anomalous dimension" of 
operator O a . A closer analysis will tell that among the operators [O a ], there must be kinetic ones for which 
the above variations would induce an extra finite "renormalization" of certain field operators, the rest must 
be characterizing interactions for which the variations would give rise to beta functions of their couplings [6~H- 
[66j|. So, Eq. (|89[) is our primitive form of renormalization group equations for a general prescription, which 
could be naturally interpreted as a "decoupling theorem" of the underlying structures that regulate the 
ultra-violet regions of our quantum field theories. Obviously, imposing different boundary conditions would 
lead to different form of d^Cidd and ^ a So a Io a ■ 

In fact, one could recast the scaling law of Eg. ((57]) into the following concise form 

{p ■ d p - d r(n) }rW(b], [m, e]; [Q]) = I^HMi K e]; [G]) 

= rW([0;p],[m,e];[Cy), (90) 

with I t Q = —md m — J2i ddCidd denoting the insertion of the full trace of the stress tensor that contains 
trace anomalies due to renormalization. This equation is just an alternative form of the Callan-Symanzik 
equation, from which some low-energy theorems of QCD follow as immediate corollaries [67). 
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V. DISCUSSIONS AND SUMMARY 



Before closing our presentation, we wish to make the following remarks about our simple strategy for 
renormalization: (1) As we introduce no artificial deformation of all the canonical structures or symmetries, 
hence not affecting spacetime dimension and Dirac algebra, thus our strategy could be applied to any held 
theory, such as supersymmetric and chiral field theories. (2) Since our ignorance about the underlying theory 
is parametrized in a general manner, the effects of regularization upon field theory and corresponding physics 
could be examined in a general manner. In particular, the physical origin of field theoretical properties like 
anomalies could be unambiguously identified and disentangled with the regularization effects in this strategy, 
surpassing the traditional interpretations that rely on effects of regularization. (3) Various operations such 
as external momenta routing and shift of loop variables are allowed in our simple strategy, sparing many 
subtleties associated these operations in the conventional approaches. Hence, one could better focus on more 
physical issues in calculations or identify true physical origins of various phenomena within field theoretical 
contexts. (4) In theoretical perspective, our approach is also useful in exploring whether or to what extent the 
underlying structures are compatible with the canonical symmetries or properties of quantum field theories, 
such as Lorentz invariance, various gauge symmetries, unitarity and so on. From a more practical viewpoint , 
our approach or similar strategies allows us to examine how our ignorance about underlying structures could 
be safely treated with the helps of these canonical symmetries. Obviously, more works need to be done both 
in the construction and applications of efficient programs using the simple strategy and in the exploration 
of various issues that are intricate with respect to the issue of regularization and renormalization. 

In summary, we demonstrated in details with QED at one-loop level how a simple strategy for renormal- 
ization could work without introducing any specific form of regularization and manipulations of ultra-violet 
infinities. Some technical operations like loop variable shifting or momenta routing that are subtle issues 
in conventional programs are shown to be of no problem. In this simple strategy, ambiguities arise instead 
of infinities, which could be further reduced by imposing appropriate Ward identities. For the Lorentz and 
gauge invariant QED at one-loop level, it is shown in a prescription-independent manner that there are 
only three ambiguities to be fixed using appropriate boundary conditions. The conventional renormalization 
programs were also analyzed from the underlying theory perspective with the rationale of manipulation of 
ultra-violet infinities being explicated. Finally, the renormalization-group-like equations arise generically as 
"decoupling theorems" of the underlying structures. 



Acknowledgement 



The project is supported in part by the National Natural Science Foundation under Grant No. 10205004 
and the Ministry of Education of China. 



Appendix A 

In this appendix, we list the convergent integrals that appear in the RHS of Eqs. (|14ll5ll6[) . 



i£(i);«*/»(p) = {-ief 



dH i [(1-0^-9^] 
P + ie 

-7" + (a ^ (3) 



(2tt)4 
-1 



7« 



A«;>,p)= {-ief 



/ + j> — m + ie 

dn 1 [(i-o9 1 -.9 kt ] 



/ + j> — m + ie 
-1 

/ + p — m + ie 

i 



(2tt)4 



-1 



l 2 + ie 

i 



"7* 



/ — m + ie 



r7 



I + j) — tn + ie ^ / + ^5 — m + ie 



in#^(g)= (-ze) 5 



(2tt)< 



Tr 



7m 
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-1 



1 — 7' 

/ + ^ — m + ie J, + f — m + ie / + g — m + ie 
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(Al) 
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-permutations of(a, /3, 7). 



(A3) 



Obviously, these diagrams result in from zero momentum insertion of photon lines into the original diagrams 
that are ill-defined, Cf. Fig. 2 for a diagrammatic representation for the case of self-energy diagram. 

As stated in Sec. II and III, other choices of routing only lead to equivalent results. For the self-energy 
diagram, we could let the external momentum p of fermion flow through the photon line, so that the RHS 
of Eq.(Ql]) reads, 



4; I 



dH 



(1 e\ {1-pT{1-pY 
\ L ~ ?J — (T^p — 



(2tt)4 



(I -p) 2 +ie 



"7k 



1 r^Tr 

/ — m + it 



(A4) 



It is nothing else but the original one with the loop momentum shifted by p. Then applying our method to 
this integration, we arrived the following, 



iE«; alt (p) 



icr 



dz [[(z(2 - 3z) + z(3z - 4)£$ + (4 + 2z(£, - l))m 
p 2 (m + zf) 



(4tt) 2 

x ln(m 2 - zp 2 ) - (1 - i)z(l - z) 



m 2 — zp 2 



+ c^i> + c„ 



(A5) 



which is equivalent to the form of the self-energy diagram as given in Eq. (|24[) after carrying out the parametric 
integrations, and the constants may differ from the those in Eq. (l24[) by a finite amount. Indeed, computing 
the difference between this routing and that given in Eq. (|2~4"l) . we have, 



(l);alt 



(4tt) 2 



Cib C-il, 



(C-l) 



Thus the two routings produce identical result after setting 

(1-0 



Cql} Cml 



(A6) 



(A7) 



which is natural from the underlying theory perspective where we should have -{(E' 1 '' 8 " — E^) = 0. Of 
course, one could also verify this point using a consistent regularization method like dimensional scheme. 
This is an illustration of the theorem of routing in case (a) given in Sec. II. 

Alternatively, we could proceed as follows: First, we could find the following from what we obtained above 
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Next note that, 
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Then from underlying theory perspective the two quantities — idpaH^ || p= o and — i9 pa S( 1 ) ;alt ||p = o should be 
identical and hence C$ = + In fact, the integral forms of these two quantities look different from 

each other, but it could be verified in any consistent regularization (e.g., dimensional regularization) that 
their difference is zero, 



dH e 2 (3-07^ 



) 4 l 2 (l 2 - m 2 ) 
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I 2 



I 2 — m 2 



= 0, 



(A13) 



therefore they are identical in any consistent regularization scheme (including the postulated underlying 
theory), and the case (a) of the theorem of routing in Sec. II is verified. 

We also note that the routing for the vertex diagram given above are so chosen that the calculations are 
less laborious, i.e., the two external fermion lines carry independent momenta so that the differentiation with 
respect to p only operates on one internal line. One could well try the following more conventional routing, 



A<, 1);a (p,p + «) = (-ief 
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(A14) 



The final results for A M should be the same after we set q = p — p in the solutions. After quite some algebras, 
we found that the vertex computed using this routing reads, 
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V + + 37m W + + 3^7m + (m 2 - 6P 2 ) 7m 

+2m(3P /1 - 2 Pfl - qj. 



(A15) 



It is immediate to see that this result exactly agrees with that given above after setting p = p + q, as asserted 
by the theorem of routing in case (b) presented in Sec. II. 

Appendix B 

Here, we consider the Ward identity for the components proportional to 1 — £ that are listed as below: 
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with 5 = m 2 /p 2 . Differentiating AS^^p; e) with respect to p^, we have, 
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Obviously, the Ward identity for the components proportional to (1 — £) is guaranteed by the following 
relation between AC^(e) and ACi(e): 



ACV(e) + AC ± (e) - 



11(^-1) 



Now, collecting all the components for C^{e) and Cj_(e): 

C„(e) = C^(e)|| ?=1 + ACV(e), C ± (c) - C ± (e)|| €=1 + AC ± (e), 

we have: 

C„(e) + Cj_(e) -- + = — — . 



In addition, we note that in Landau gauge, 
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as it is well known that the vertex is ultra-violet definite in this gauge. 
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Appendix C 



The differentiation of other diagrams with respect to pi- a yields, 
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with 
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r 6 , 




(C2) 
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with the same convention of notations. Here we have deliberately denote the integration constants in the 
way that seems diagram-dependent, i.e., C ; 4 7 may differ among the six diagrams. The gauge invariance will 
require that they must be equal to each other. One could also arrive at the same conclusion by noting that 
the loop integral for each diagram is identical after setting all the external momenta to zero, for example, 



Thus the solutions: 
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Appendix D 



Now we demonstrate the integration by parts on the logarithmic term InAi: 
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with 



Ai = Ai|| 2=a , D lz = - I dx I dy I dz 



Ai 

~ r 1 , r, y 2 a„Ax - 1 r\ x 3 (i-2x) 

D.^-^dxJ^^-, D lx = -J a dx- r — T - m , 

D 1 =D lx + D ly + D lz . (D2) 

Here Di is definite and hence absorbed into lj in Eq. (|4"3")l . that is, f 1 ! = r x + D\, Then we arrive at the 
form of Ti-^pa given in Eq. (|43[) . In the same fashion, we have, 

r 1 r x f v ~ 1 m 2 

/ dx / dy / ^lnA fc = -In- + D k , k = 2,3,4,5,6, (D3) 

JO JO JO ° ( - y fc;47 

with [Dfe, fc = 2, 3, 4, 5, 6] being definite functions of external momenta. 
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FIG. 1: The elementary one-loop diagrams in QED 
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FIG. 2: The graphical representation of Eq.(14). 



